Conduction channels of superconducting quantum point contacts 
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Atomic quantum point contacts accommodate a small number of conduction channels. Their number 
N and transmission coefficients {Tn} can be determined by analyzing the subgap structure due to 
multiple Andreev reflections in the current-voltage {IV) characteristics in the superconducting state. 
With the help of mechanically controllable break-junctions we have produced Al contacts consisting 
of a small number of atoms. In the smallest stable contacts, usually three channels contribute to the 
transport. We show here that the channel ensemble {Tn} of few atom contacts remains unchanged 
up to temperatures and magnetic fields approaching the critical temperature and the critical field, 
respectively, giving experimental evidence for the prediction that the conduction channels are the 
same in the normal and in the superconducting state. 



I. INTRODUCTION 

An atomic size contact between two metallic electrodes 
can accommodate only a small number of conduction 
channels. The contact is thus fully described by a set 
{Tn} = {Ti,T2, ...Tn} of transmission coefficients which 
depends both on the chemical properties of the atoms 
forming the contact and on their geometrical arrange- 
ment. Experimentally, contacts consisting of even a sin- 
gle atom have been obtained using both scanning tunnel- 
ing microscope and break-junction techniques The 

N 

total transmission D~ ^ Tn of a contact is deduced 

n=l 

from its conductance G measured in the normal state, us- 
ing the Landauer formula G ~ GqD where Go = 2e^/ft- 
is the conductance quantum 

Experiments on a large ensemble of metallic contacts 
have demonstrated the statistical tendency of atomic-size 
contacts to adopt configurations leading to some pre- 
ferred values of conductance. The actual preferred values 
depend on the metal and on the experimental conditions. 
However, for many metals, and in particular 'simple' ones 
(like Na, Au...) which in bulk are good 'free electrons' 
metals, the smallest contacts have a conductance G close 
to Go [0J|^ • Statistical examinations of Al point con- 
tacts at low temperatures yield preferred values of con- 
ductance at G = 0.8 Go, 1.9 Go, 3.2 Go and 4.5 Go §, 
indicating that single-atom contacts of Al have a typi- 
cal conductance slightly below the conductance quantum. 
Does this mean that the single-atom contacts correspond 
to a single, highly transmitted channel (T = 0.8)? This 
question cannot be answered solely by conductance mea- 
surements which provide no information about the num- 
ber or transmissions of the individual channels. 

However, it has been shown that the full set {Tn} is 
amenable to measurement in the case of superconducting 
materials j7| by quantitative comparison of the measured 
current-voltage {IV characteristics) with the theory of 



multiple Andreev reflection (MAR) for a single channel 
BCS superconducting contact with arbitrary transmis- 
sion T, developed by several groups for zero tempera- 
ture = and zero magnetic field H — ||^-[ll|]. Al- 
though the typical conductance of single-atom contacts 
of Al (G ~ 0.8 Go) is smaller than the maximum possible 
conductance for one channel, three channels with trans- 
missions such that Ti + T2 + ~ 0.8 have been found 

Moreover, there exist other physical properties which 
are not linear with respect to {Tn} as e.g. shot noise 
1 16 1, conductance fluctuations jl^ , and thermopower 
1 18 1, which also give information about the {Tn} of a 
contact. Although it is not possible to determine the 
full set of transmission coefficients with these properties, 
certain moments of the distribution and in particular the 
presence or absence of partially open channels can be 
detected. Recent experiments have shown that normal 
atomic contacts of Al with conductance close to Go con- 
tain incompletely open channels p^ ] in agreement with 
the findings in the superconducting state [[7|. 

In previous work we have shown ||l^ , |l^ how the con- 
duction channels of metallic contacts can be constructed 
from the valence orbitals of the material under investiga- 
tion. In the case of single-atom contacts the channels are 
determined by the valence orbitals of the central atom 
and its local environment. In particular for Al the chan- 
nels arise from the contributions of the s and p valence 
bands. To the best of our knowledge it has never been 
observed in contacts of multivalent metals that a single 
channel arrives at its saturation value of T = 1 before 
at least a second one had opened. Single-atom contacts 
of the monovalent metal Au transmit one single channel 
with a transmission < T < 1 depending on the partic- 
ular realization of the contact 



14 



151 



From the theoretical point of view no difference be- 
tween the normal and superconducting states is expected, 
because {i) according to the BCS theory |l^] the elec- 
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tronic wave functions themselves are not altered when 
entering the superconducting state, but only their occu- 
pation and (a) MAR preserves electron-hole symmetry 
and therefore does not mix channels ]20| , pT[ |. 

Experimental evidence for the equivalence of the nor- 
mal and superconducting channels can be gained by trac- 
ing the evolution of the IV curves from the supercon- 
ducting to the normal state in an external magnetic field 
and/or higher temperatures and comparing them to the 
recent calculations by Cuevas et al. of MAR in single 
channel contacts at finite temperatures and includ- 
ing pair breaking due to magnetic impurities or magnetic 
field ||. 

We show here that the channel ensemble {T^} of few 
atom contacts remains unchanged when suppressing the 
superconducting transport properties gradually by rais- 
ing the temperature or the magnetic field up to temper- 
atures and magnetic fields approaching the critical tem- 
perature and the critical field, respectively. Although it is 
not possible to measure the full channel ensemble above 
the critical temperature or field, respectively, no abrupt 
change is to be expected since the phase transition (as 
a function of temperature) is of second order. Because 
the determination of the channel ensemble relies on the 
quantitative agreement between the theory and the ex- 
perimental IVs, we concentrate here on the case of Al 
point contacts since we expect for this material the BCS 
theory to fully apply. 
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FIG. 1. Measured IV curves (symbols) of four different 
atomic contacts with G ~ 0.9 Go at < 50 mK and best 
numerical fits (lines). The {Tn} and total transmissions D 
obtained from the fits are: (a) {0.900,0.108}, D = 1.008; 
(b) {0.802,0.074}, D = 0.876; (c) {0.747,0.168,0.036}, D = 
0.951; (d) {0.519, 0.253, 0.106}, D = 0.878. Current and volt- 
age are in reduced units, the current axis is normalized with 
respect to the total conductance measured by the slope of the 
IV at high voltages eV > 5A. Not all measured data points 
are shown. The measured gap was A/e — (184 ± 2) /^V. 
Left inset: Theoretical IVs for a single channel superconduct- 
ing contact for different values of T (from bottom to top: 0.1, 
0.6, 0.9, 0.99, 1) after 0. Right inset: typical total transmis- 
sion traces measured as a function of the electrode distance 
d at eV > 5 A, while opening the contact at around 6 pm/s. 
The d axis has been calibrated by the exponential behavior 
in the tunnel regime. 



II. TRANSPORT THROUGH A 
SUPERCONDUCTING QUANTUM POINT 
CONTACT 

The upper left inset of Fig. ^ shows the theoretical IVs 
by Cuevas et al. |10| for zero temperature G = and zero 
field H = 0. A precise determination of the channel con- 
tent of any superconducting contact is possible making 
use of the fact that the total current I{V) results from 
the contributions of N independent channels: 



N 



/(y) = ^z(v,T„ 



(1) 



This equation is valid as long as the scattering matrix 
whose eigenvalues are given by the transmission coeffi- 
cients is unitary, i.e. the scattering is time independent. 
The i{V, T) curves present a series of sharp current steps 
at voltage values V = 2/S./me, where m is a positive inte- 
ger and A is the superconducting gap. Each one of these 
steps corresponds to a different microscopic process of 
charge transfer setting in. For example, the well-known 
non-linearity at eV = 2A arises when one electronic 
charge (m = 1) is transferred thus creating two quasi- 
particles. The energy eV delivered by the voltage source 
must be larger than the energy 2A needed to create 



the two excitations. The common phenomenon behind 
the other steps is multiple Andreev reflection (MAR) 
of quasiparticles between the two superconducting reser- 
The order m = 2 , 3 , . . . , of a step corresponds 



25 1 



to the number of electronic charges transferred in the un- 
derlying MAR process. Energy conservation imposes the 
threshold meV > 2A for each process. For low trans- 
mission, the contribution to the current arising from the 
process of order m scales as T™. The contributions of all 
processes sum up to the so-called "excess current" the 
value of which can be determined by extrapolating the 
linear part of the IVs well above the gap v > 5A down 
to zero voltage. As the transmission of the channel rises 
from to 1, the higher order processes grow stronger 
and the current increases progressively. The ensemble of 
steps is called " subharmonic gap structure" , which was 
in fact discovered experimentally ]2^ , has been exten- 
sively studied in superconducting weak links and tunnel 
junctions with a very large number of channels [E7],E8|. 



III. EXPERIMENTAL TECHNIQUES 

In order to infer {Tn} from the IVs, very stable atomic- 
size contacts are required. For this purpose we have 
used micro-fabricated mechanically controllable break- 
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FIG. 2. Three point bending mechanism. The distance be- 
tween the two counter-supports is 12 mm, and the substrate 
is 0.3 mm thick. The micrograph shows a suspended Al mi- 
crobridge. The insulating polyimide layer was etched to free 
the bridge from the substrate. The third panel displays the 
wiring of the experimental setup. 



junctions Our samples are 2/xm long, 200 nm thick 
suspended microbridges, with a 100 nm x 100 nm con- 
striction in the middle (cf. Fig. The bridge is broken 
at the constriction by controlled bending of the elastic 
substrate mounted on a three-point bending mechanism. 
A differential screw (100 /xm pitch) driven by a dc-motor 
through a series of reduction gear boxes, controls the mo- 
tion of the pushing rod that bends the substrate (Fig. H). 

The geometry of the bending mechanism is such that 
a 1 /im displacement of the rod results in a relative mo- 
tion of the two anchor points of the bridge of around 
0.2 nm. This was verified using the exponential depen- 
dence of the conductance on the interelectrode distance 
in the tunnel regime. This very strong dependence was 
used to calibrate the distance axis to an accuracy of about 
20 %. The bending mechanism is anchored to the mixing 
chamber of a dilution refrigerator within a metallic box 
shielding microwave frequencies. The bridges are broken 
at low temperature and under cryogenic vacuum to avoid 
contamination. 

The IV characteristics are measured by voltage biasing 
with U — Udc the sample in series with a calibrated resis- 
tor Rs = 102.6 kfl and measuring the voltage drop across 
the sample (giving the V signal) and the voltage drop 
Vs = IRs across Rg (giving the / signal) via two low-noise 
differential preamplifiers. The differential conductance is 
measured by biasing with U = Udc + Ui cos{2tt ft) using a 
lock-in technique at low frequency / < 200 Hz. All lines 
connecting the sample to the room temperature electron- 
ics are carefully filtered at microwaves frequencies by a 
combination of lossy shielded cables |Q , and microfab- 
ricated cryogenic filters [0. The cryostat is equipped 
with a superconducting solenoid allowing to control the 
field (IqH at the position of the sample within 0.05 mT. 
After having applied a magnetic field and before taking 
new H = data we demagnetize carefully the solenoid. 
The temperature is monitored by a calibrated resistance 
thermometer thermally anchored to the shielding box. 
The absolute accuracy of the temperature measurement 
is about 5%. 



IV. DETERMINATION OF THE CHANNEL 
TRANSMISSIONS 

Pushing on the substrate leads to a controlled open- 
ing of the contact, while the sample is maintained at 
Q < 100 mK. As found in previous experiments at higher 
temperatures, the conductance G decreases in steps of 



the order of Gg , their exact sequence changing from open- 
ing to opening (see right inset of Fig. |l|) . The last conduc- 
tance value before the contact breaks is usually between 
0.5 and 1.5 Gq. 

Figure |l| shows four examples of IVs obtained at 
O < 50 mK on last plateaux of two different Al samples 
just before breaking the contact and entering the tun- 
nel regime. The curves differ markedly eventhough they 
correspond to contacts having the same conductance of 
about G c^i 0.9 Go within 10%. The existence of IVs 
with the same conductance but different subgap struc- 
ture implies the presence of more than one channel with- 
out further analysis. In particular, the examples shown 
here demonstrate that although they would correspond 
to the first maximum in the conductance histogram [||, 
they do not transmit a single channel, but at least two, 
with a variety of transmissions. 

In Fig. |l| we also show the best least-square fits ob- 
tained using the numerical results of the Q = theory 
of Cuevas et al. (l^ . The fitting procedure decomposes 
the total current into the contributions of eight indepen- 
dent channels. Channels found with transmissions lower 
than 1% of the total transmission were neglected. When 

< 3, this fitting procedure allows the determination of 
each Tn with an accuracy of 1% of the total transmis- 
sion D. For contacts containing more channels only the 
2 or 3 dominant channels (depending on their absolute 
value) can be extracted with that accuracy. Details of 
the fitting procedure are published in ||^. 



V. IVS OF AL POINT CONTACTS AT HIGHER 
TEMPERATURES 

The right inset of Fig. ^ displays the evolution of the 
IV of contact (a) from Fig. |l| for three different temper- 
atures below the critical temperature 8c — 1.21 K (each 
trace is offset by 0.5 for clarity). When the temperature 
is increased the subgap structure is slightly smeared out 
due to the thermal activation of quasiparticles, and the 
position of the current steps is shifted to smaller volt- 
ages due to the reduction of the superconducting gap. 
Although the IVs are very smooth due to the dominance 
of an almost perfectly open channel, up to eight MAR 
processes are distinguishable in the dl/dV. The solid 
lines are calculated with the same set of transmissions 
{Ti = 0.900, T2 = 0.108} for the temperatures indicated 
in the caption using the BCS dependence of the super- 
conducting gap and the Fermi function for the respective 
temperature |22 2^. The quality of the fit does not vary 
with temperature. 

In order to further interpret the data, we plot in Fig. ^ 
the theoretical zero temperature single channel dl/dVs 
for the same transmissions as in the inset in Fig. |^ in 
two different manners. In the left panel they are plotted 
as a function of eV/ A showing the different shapes and 
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FIG. 3. Differential conductance dl/dV as a function of 
voltage measured at = 47 mK (bottom), 640 mK (mid- 
dle), and 810 mK (top) for a contact with total conductance 
D — 1.008. The two upper curves are offset by 2 each. 
The solid lines are fits using the same transmission ensemble 
{Ti = 0.900, T2 = 0.108} for all three temperatures. Right in- 
set: corresponding IVs offset by 0.5 each. Left inset: position 
of the m = 1 maximum in the d//d1^ curves. Data has been 
taken on three different contacts with G varying from 1 to 2 
Go. The line is the temperature dependence of the supercon- 
ducting gap according to BCS theory for A{0 = 0) = 184/iV. 
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FIG. 4. Calculated dl/dV characteristics (normal- 
ized to the conductance) of a single channel quantum 
point contact for transmission 0.1 (solid), 0.6 (dashed), 0.9 
(dash-dotted), 1.0 (dotted) as a function of eV/ A (left panel) 
and m' = 2A/eV (right panel). The maximum of the curve 
for T = 0.1 arrives at 12. 6G. The amplitude, the shape 
and the exact position of the individual maxima depend on 
the transmission. Only for small T the maxima are close to 
V = 2 A/me, for higher transmission they are shifted to lower 
voltages. The subgap structure is most pronounced for inter- 
mediate T. 

with a slope corresponding to D within 1%. 



amplitudes of the individual MAR processes for varying 
transmission. A small T gives rise to narrow conduc- 
tance spikes, whereas a higher T yields round maxima 
at voltages V < 2 A/me and pronounced minima close 
to the sub- multiple values V = 2 A/me. This behavior 
is clearly visible in the right panel where the differential 
conductance is plotted as a function of the generalized 
order of the MAR process m' = 2A/eV. For small T 
the onset of the MAR processes is equidistant with spac- 
ing 2A/eV and their amplitudes decrease very rapidly 
with m'. For high transmission the position of the max- 
ima is progressively shifted to higher m' values, while 
the minima correspond approximately to integer values 
of m'. The experimental data of Fig. || display a mixed 
character of high T and low T behavior, because of the 
presence of the two extreme channels with Ti = 0.90 and 
T2 = 0.108. The value of the temperature dependent 
gap A(8) can be determined by the peak of the m = 1 
process, while the rest of the dl/dV is dominated by the 
widely open channel Ti. In the left inset of Fig. || we 
plot the position of the m = 1 maximum as a function 
of the temperature. Also shown are data taken on dif- 
ferent contact configurations of the same sample. The 
development of the peak position follows the BCS gap 
function Abcs(0) which is plotted as a solid line in the 
same graph. We have verified for contacts with different 
conductances ranging from the tunnel regime to several 
Go that the IVs can be described by the same channel 
distribution (with restricted accuracy due to less pro- 
nounced MAR features) up to the critical temperature. 
When exceeding 8c the IV characteristics become linear 



VI. IVS OF AL POINT CONTACTS WITH 
EXTERNAL MAGNETIC FIELD 

Fig. H shows the evolution of the subgap structure with 
applied magnetic field. The traces are offset for clarity. 
In our experiment the field is applied perpendicular to the 
film plane. As the field size is increased, the excess cur- 
rent is suppressed, the current steps are strongly rounded 
and the peak positions are shifted to lower voltages. For 
fields larger than w 5.0 mT no clear sub-multiple current 
steps are observable. When a field of ^qHc = 10.2 mT 
close to the bulk critical field of Al yUoi?c,buik = 9.9 mT is 
reached the IV becomes again linear with a slope corre- 
sponding to the sum D of the transmissions determined 
in zero field. When lowering the field again to H = 
we recover the same subgap structure as before. When 
reversing the field direction the same IV is observed for 
the same absolute value of the field, which proves that 
there is no residual field along the field axis. Effects of 
the earth magnetic field or spurious fields in different di- 
rections however cannot be excluded. 

An external magnetic field suppresses superconductiv- 
ity because it acts as an effective pair breaking mecha- 
nism . Since a magnetic field breaks the electron-hole 
symmetry, the IVs due to MAR are modified as demon- 
strated by Zaitsev and Averin p^ . Strictly speaking, the 
conduction channels could be different with and without 
magnetic field. 

A quantitative description of the influence of the mag- 
netic field is difficult because of the complicated shape of 
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FIG. 5. Measured (symbols) and calculated (lines) IVs ol 
trace (d) of Fig. ^ with transmission set {0.519, 0.253, 0.106} 
for magnetic fields of fioH — 0, 2.1, 5.4, 5.9, 7.7, 9.1 and 
10.9 mT (from bottom to top). The curves with magnetic 
field are offset by 0.2 each. Left inset: pair amplitude A 
(dotted), spectral gap Q (solid) (according to ^^) and posi- 
tion of the m — 1 maximum of the corresponding calculated 
d//dVs (symbols) in units of the gap parameter at zero field 
as a function of the pair breaking parameter F. Right inset: 
Pair breaking parameter F used for calculating the IVs and 
dl/dVs as a function of applied field (symbols) and prediction 
according to (hne). The IV a.t H = was measured 

at the beginning and the end of the series of measurements 
and the {Tn} was found to be reproducible. The critical field 
at the measuring temperature of = 25 mK was determined 
to be ^ioHc = 10.2 mT. 



the samples. The point contact spectra are sensitive to 
the superconducting properties at the constriction. Since 
the pair breaking parameter F — ?i/(2TpbA) (rpb is the 
pair-breaking time) due to an external magnetic field is 
geometry dependent |32| a complete description needs 
to take into account the exact shape of the sample on 
the length scale of the coherence length ^. We estimate 
for our Al films in the dirty limit ^ = {hDei/2Ay^^ = 
280 nm where Dc\ = vpl/3 — 0.042 m^/s is the electronic 
diffusion constant. 

Due to the finite elastic mean free path / « 65 nm 
(determined by the residual resistivity ratio RRR = 
i?(300K)/i?(4.2K) « 4) of the evaporated thin film the 
penetration depth is enhanced and it is comparable to the 
sample thickness. The fact that all signature of supercon- 
ductivity is destroyed at the bulk critical field indicates 
that the geometry of the sample does not play a domi- 
nant role, but that ^ is the most important length scale. 
We therefore describe the influence of the magnetic field 
along the lines of Skalski et al . using a homogeneous 
F given by the expression [3^j35| ] 



F = 



6hA 



(2) 



where the effective width of the film w = 280 nm ~ ^ 
is limited by the coherence length. Superconductivity is 
completely suppressed when F = 0.5. In order to ob- 
tain the IV curves for one channel in an external mag- 
netic field, the BCS density of states in the theory of Ref. 
pOf is replaced by the corresponding expressions given in 
Refs. which include the effect of a pair breaking 

mechanism. 

Contrary to the influence of higher temperatures, the 
magnetic field rounds the density of states and the An- 
dreev reflection amplitude. The rounding is a conse- 
quence of the fact that the pair amplitude A and the 
spectral gap fl of the density of states (i.e. the energy 
up to which the density of states is zero) differ from each 
other when time reversal symmetry is lifted. The posi- 
tion of the m = 1 maximum of the d//dT^ does not give 
an accurate estimation of F and it is necessary to fit the 
whole IV. In the left inset we display the evolution (as a 
function of F) of A, and the position of the maximum 
conductance of the m = 1 process of the contact whose 
IVs are shown in Fig. |[ The functional dependence of 
the latter is not universal but depends on the distribution 
of transmissions. 

It turns out that the structure in the experimental data 
is more rounded than in the calculated curves indicating 
the limitations of the model used here. A reasonable 
agreement between the experimental data and the model 
is found when determining F such that the excess cur- 
rent is correctly described. The solid lines in Fig. |^ are 
calculated with the transmission ensemble determined at 
zero field for values of F given in the inset. These values 
correspond nicely to the predicted quadratic behavior of 
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Eq. g. It was not possible to achieve better agreement 
between the measured and calculated IVs by altering the 
channel ensemble, supporting again that the conduction 
channels are not affected by superconductivity, nor by its 
suppression p6[ |. 

We have demonstrated here, that it is possible to drive 
a particular contact reproducibly into the normal state 
and back to the superconducting state without changing 
{Tn}. We stress the high stability of the setup necessary 
for maintaining a particular contact stable during the 
measurement series. 



VII. CONCLUSIONS 

We have reported measurements and the analysis of 
multiple Andreev reflection in superconducting atomic 
contacts demonstrating that the conduction channel en- 
semble of the smallest point contacts between Al elec- 
trodes consists of at least two, more often three channels. 
We have verified that the channel ensemble remains un- 
changed when suppressing the superconductivity gradu- 
ally by increasing the temperature or applying a mag- 
netic field. This result strongly supports the expected 
equivalence of conduction channels in the normal and in 
the superconducting state and agrees with the quantum 
chemical picture of conduction channels. The latter sug- 
gests that the conduction channels are determined by the 
band structure of the metal and therefore their transmis- 
sions vary significantly only on the scale of several eV. 
Superconductivity, which opens a spectral gap for quasi- 
particles of the order of only ~ meV does not modify the 
channels and is therefore a useful tool to study them. 
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